Direct measurement of the ^^N(p,7)^^0 5"- factor 
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We have measured the ^'*N(p, 7)^^0 excitation function for energies in the range E^"* = 155- 
524 keV. Fits of these data using _R-matrix theory yield a value for the S-factor at zero energy of 
1.64 ± 0.07 (stat) ±0.15 (sys) keV-b, which is significantly smaller than the result of a previous 
direct measurement. The corresponding reduction in the stellar reaction rate for ^*N(p, 7)^^0 has 
a number of interesting consequences, including an impact on estimates for the age of the Galaxy 
derived from globular clusters. 
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Stars produce energy from the conversion of hydrogen 
into hehum primarily by the p-p chains and by the CN- 
cycle. The latter is the dominant energy source for stars 
somewhat more massive than the sun, but all stars will 
produce energy via the CN cycle at the end of their main- 
sequence lifetimes, and while on the red-giant branch. 
At the burning temperatures characteristic of these evo- 
lutionary stages (T « 0.02-0.055 GK), the i'^N(p, 7)1^0 
reaction is the slowest CN reaction and thus it regulates 
the rate of nuclear energy generation. The power liber- 
ated by the CN cycle and the amount of helium produced 
has a direct connection to the luminosity observed at the 
transition between the main sequence and the red-giant 
branch, and on the luminosity of the horizontal branch. 
Both of these quantities play a role in determining the 
ages of globular clusters [1, 2]. Also, since it helps to 
constrain the temperature and density profiles in the H- 
burning shell, ^'^N{p, 7)^^0 will affect nucleosynthesis be- 
yond the CN cycle during the red-giant stage. 

The thermonuclear reaction rate can be obtained from 
the astrophysical S-factor, defined as 



(1) 



where Ecm is the energy in the center of mass, (T(Ecm) is 
the reaction cross section, and 77 is the Sommerfeld pa- 
rameter. The accepted S-factor for the ^''N(p, 7)^^© re- 
action is based on the measurements of Schroder et al. [3] 
(hereafter Sch87) who quoted a zero-energy 5-factor of 
5(0) = 3.20 ± 0.54 keV-b. About half of this S-factor re- 
sults from transitions into the bound state at Ex = 6.79 
MeV (J''=3/2+) in ^^O. Most of the remainder arises 
from capture to the ground state (J'^=l/2^), including 
a significant contribution from the tail of the 6.79-MeV 
state (corresponding to a subthreshold state at Ecm = 
-0.504 MeV). However, uncertainty about the inferred 
width of the subthreshold state led to a recommendation 
of S[Q) — 1.5-4.5 keV-b [4], and a subsequent re-analysis 
of the data of Sch87 by Angulo and Descouvemont [5] re- 
ported S(0) = 1.77 ± 0.20 keV-b. A measurement of the 



lifetime of the 6.79-MeV state [6] strongly favored the 
lower range of values for S{Q). More recently, another 
re-analysis of the Sch87 data using measured asymp- 
totic normalization coefficients (ANCs) as constraints 
produced a consistent resuh, 5(0) 1.70±0.22 keV-b [7]. 
It should be noted that in the vicinity of the 0.259- and 
0.985-MeV resonances, Sch87 present yield data rather 
than cross sections and that errors were made in their cor- 
rections for coincident summing [8]. Thus, a number of 
these points should not have been included in the analy- 
ses described above. Nonetheless, evidence favors a value 
for the 5-factor that is 20-40% smaller than previously 
thought for the temperatures of interest. Because of the 
importance of the ^'*N(p, 7)^^© reaction in determining 
the power liberated by the CN cycle, even a change of this 
magnitude impacts all of the issues mentioned above. In 
view of the importance of this reaction for stellar astro- 
physics, we have carried out a new measurement designed 
to more accurately determine the low-energy 5-factor. 

We measured an excitation function for the 
^'*N(p, 7)^^0 reaction at the Laboratory for Exper- 
imental Nuclear Astrophysics (LENA), located at the 
Triangle Universities Nuclear Laboratory. A 1-MV Van 
de Graaff accelerator provided proton beams at labora- 
tory energies between 155 and 524 keV, and with beam 
currents of 100 - 150 /lA. The beam entered the target 
chamber through a copper tube, extending to less than 1 
cm from the surface of the target, which was cooled using 
chilled, de-ionized water. The copper tube was cooled by 
a LN2 reservoir to trap potential target contaminants. 
In order to suppress the emission of secondary electrons 
from the target, permanent magnets were located at 
the end of the tube along with an electrode biased to 
-300 V . Targets were fabricated by implanting nitrogen 
ions into 0.5 mm-thick tantalum backings. The backings 
were prepared for implantation by etching in an acid 
solution [9] to remove surface impurities that could give 
rise to beam-induced background. Implantation energies 
of 20, 40 and 110 keV were used to produce targets that 
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were 5, 10 and 18-kcV thick, respectively, as measured 
at the Ecm = 0.259-MeV resonance [10]. TypicaUy, 
their composition and thickness remained stable over 
accumulated doses of 20-25 C. The condition of the 
target was checked periodically by measuring a yield 
curve for the 0.259-MeV resonance. The stoichiometry 
(Ta/N = 0.718±0.076) was measured via Rutherford 
backscattering and was consistent with published val- 
ues [11]. Prom 32 independent measurements of the 
thick-target yield, we obtain lo^ = 0.0135±0.0012 eV for 
the 0.259-MeV resonance, which is in good agreement 
with the previous value of 0.014±0.001 eV [12]. Our 
uncertainty is purely systematic and arises primarily 
from the uncertainty that we assign to the stopping 
powers (calculated using SRIM2000 [13]). By compari- 
son, the statistical uncertainty (0.4%) is negligible. We 
have also measured the branching ratios for the decay 
of the 0.259-MeV resonance (Usted in Table I). This 
measurement was carried out with the HPGe detector 
moved to a distance of 23 cm from the target, which 
reduced coincidence summing effects to a negligible level. 
It should be noted that our value for the ground-state 
transition differs significantly from the value of Sch87. 

Gamma rays were detected using a 135% HPGe de- 
tector placed at 9iab = 0° and at a distance of 9 mm 
from the target. The energy calibration and absolute 
photopeak efficiency were established using radioactive 
sources and the decays from well-known resonances in 
the i4n(p, 7)1^0, 26Mg(p,7)27Al and ^'^ M{p,-ff^Si re- 
actions. The total efficiency (needed for summing cor- 
rections) was calculated using MCNP [14] and normal- 
ized to source data. A 35.6-cm dia. x 40.6-cm long an- 
nulus of Nal(Tl) enclosed both the target and Ge crys- 
tal. This detector geometry allowed us to record three 
types of events: Ge singles, Ge-Nal coincidences, and Ge 
signals without a corresponding event in the Nal detec- 
tor within 5 /is. In the latter mode, the Nal served as 
a cosmic-ray veto while also suppressing events arising 
from 7-ray cascades, which proved useful for detecting 
the ground-state transition. Because of this fixed geom- 
etry, we could not measure angular distributions for the 
primary transitions. However, for our energy range, the 
primary transitions to the groimd, 5.18 and 6.18-MeV 
states were calculated to be isotropic to better than 1%. 
In addition, all secondary transitions were calculated to 
be isotropic. The only transition expected to have a sig- 
nificant angular distribution was to the 6.79-MeV state, 
which necessitated the use of the (isotropic) secondary 
transition alone to determine the cross section. The ra- 
tio of primary-to-secondary yields was consistent with 
the expected angular- distribution coefficient for an El 
transition, a2=-l. Above the 0.259-MeV resonance, the 
secondary yields had to be corrected for the possibility 
that the incident proton could lose enough energy in the 
implanted target layer to undergo resonant capture into 
this state, which would produce some of the same sec- 



ondaries. The signature of this effect was the presence 
of 7-rays from the de-excitation of the resonance, which 
also provided the means for removing their contribution 
to the secondaries. Because we used thin targets above 
the resonance, these corrections were small (7% on aver- 
age for the dominant transition to the 6.79-MeV state). 



TABLE I: Brandling; ratios for the dcc:ay of the 0.2r)9-MeV 



resonance 



transition 


this work 


Sch87 


ground state 


0.0170 ± 0.0007 


0.035 ±0.005 


5.18 MeV 


0.173 ±0.002 


0.158 ±0.006 


6.18 MeV 


0.583 ± 0.005 


0.575 ± 0.004 


6.79 MeV 


0.227 ± 0.003 


0.232 ± 0.006 



The effective center-of-mass energies, Eeff, were de- 
rived from the energies of the primary 7-rays (corrected 
for Doppler shifts) as well as from the beam energies 
and measured target thicknesses. The second procedure 
was an iterative calculation, which first assumed a lin- 
ear S-factor over the width of the target and then used 
the resulting S-factor to recalculate the result. These 
two techniques were found to yield consistent values. We 
also checked this procedure by measuring the S-factor 
with targets of different thicknesses, but at the same pre- 
dicted values for Eeff. These measiircmcnts prodiiccd S- 
factors that agreed to within experimental errors. The 
7-ray yields at each energy were corrected for coincident 
summing, which was a significant effect. In fact, for ener- 
gies below Eeff = 187 keV, the measured intensity at the 
energy of the ground state transition was entirely consis- 
tent with summing of cascades involving the higher-lying 
states. For the 6.79-MeV transition, the correction for 
summing-out was about 22% These calculations included 
the (calculated) angular correlation between the primary 
and secondary 7-rays, which was a comparatively small 
effect since the product of the attenuation coefficients 
for the primary and secondary transitions amounted to 
about 0.15. 

The three strongest transitions, in order of impor- 
tance, arc those to the 6.79-McV state, the ground state 
(g.s.), and the 6.18-MeV state (J''=3/2-). These have 
been analyzed in the framework of the J?-matrix model 
[15], including the external contributions to radiative 
capture; details can be found in Ref. [5]. These cal- 
culations have been performed with ii-matrix radii of 
a = 3.5,4,4.5,5,5.5, and 6 fm. To fit the 6.79-McV 
transition, we have included the 0.259-MeV resonance 
[J'^ = 1/2+) and the non-resonant contributions for 
J'^ = l/2~,3/2~ and 5/2~ channel spins, calculated 
from the ANC, C, of the final state. The proton width 
Tp of the resonance was left as free parameter, while its 
gamma width was fixed for all fits to F-^ = 9.2 ±0.1 meV. 
This value was obtained from our measured resonance 
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FIG. 1: _R-matrix fits of the present capture data to the 6.79 

MeV state (closed circles). The solid lino shows the fit for 
a = 5 frn; the lower (upper) dashed lines are for a = 4.5(5.5) 
fin. The corrected data from Sch87 [8] are also shown (open 
squares) but are not used in the fit. 

strength and the relevant branching ratio from Table I. 
The ANC was also left as a free parameter. The results 

of the fits arc given in Tabic II. Figure 1 shows the fits to- 
gether with the present experimental data (the notation 
RC denotes "radiative capture"). The error bars reflect 
the statistical uncertainty for each point. The data for 
this transition from Sch87 (corrected for summing effects 
[8]) are also shown in the figure for comparison, but were 
not used in the fits. As can be seen in the figure, our 
new measurement extends the S-factor to lower energies, 
with higher accuracy. 



TABLE II: Results of the ii-matrix fits for the 6.79 MeV 
transition 



a 


C 


7^ 


Tp 


SiO) 




(fm) 


(fm-V2) 


(MeV) 


(keV) 


(keV-b) 


{N = 29) 


3.5 


4.4 


1.9 


0.97 


1.1 


1.4 


4 


4.5 


1.2 


0.96 


1.1 


1.4 


4.5 


4.5 


0.8 


0.94 


1.1 


1.5 


5 


4.6 


0.5 


0.93 


1.1 


1.5 


5.5 


4.6 


0.4 


0.91 


1.1 


1.6 


6 


4.7 


0.3 


0.89 


1.2 


1.7 



The fits are weakly sensitive to the i?-matrix radius a, 
as expected. The resulting values for C and Fp are also 
not very sensitive to a. On the other hand, for a given 
a, the surface is relatively flat, indicating that several 
sets of {C,Tp) parameters are possible, but the variations 
are within the adopted uncertainty. We recommend C = 
4.5 ± 0.1 fm-i/2 and Tp = 0.94 ± 0.03 keV. The former 
value is in good agreement with an averaged value of the 
results from Mukhamcdzanov et al. [7] and Bertone et 
al. [16] (C = 4.7 ± 1.9 fm^^/^), while Tp is in very good 
agreement with the results of Sch87. For the S-factor at 
zero energy, we have obtained S{0) = 1.10 ± 0.05 keV-b. 
The uncertainties were calculated from the i?-matrix fits 



for a = 3.5 — 6 fm that have — Xmin + 1- 

The i?-matrix analysis of the g.s. transition is 
complicated by the presence of several contributing 
states. These involve the El contributions from the 
1/2+ resonance at 0.259 MeV, from the two 3/2+ 
resonances at 0.985 and 2.187 MeV, and from the 3/2+ 
subthreshold state at Ecm = —0.504 MeV. In order to 
account for higher-energy resonances, we have included 
a background pole {£ = 1) at 4 MeV with J'^=3/2+ 
and with rp=4 MeV. The location and width of this 
state do not change S(0) by more than 15%. Because 
of this large number of parameters, it was necessary to 
reduce the number of free parameters in the fit. The 
proton and gamma widths for the 0.985- and 2. 187- MeV 
resonances were taken from Ref. [5]. We have also used 
= 0.69 ± 0.03 mcV for the 0.259-MeV resonance, 
obtained from our value for 0^7 and the branching ratio 
to the ground state. The ANC is taken from Table II 
(for each value of a), while the reduced width of the 
subthreshold state 7^ is obtained from Eq. (14) of 
Ref. [5] and is also given in Table II. The number of 
free parameters is thus reduced to three: for the 
subthreshold state, Fp for the 0.259-MeV resonance, and 
the gamma width of the background pole, F^^^. The 
results of the fits are given in Table III. The fit for a = 5 
fm is shown in Figure 2. For comparison, the data of 
Sch87, corrected for summing effects [8] , are also shown, 
but they were not used in the fits. Note that these 
points are now restricted to Ecm > 324 keV whereas 
the present data extend down to Ecm = 187 keV. The 
proton width of the 0.259-MeV resonance, Fp = 1.3 ±0.1 
keV, is larger than the values of Table II obtained from 
the fits of the 6.79 MeV transition, but the error is 
larger here because of the more uncertain value for the 
branching ratio and the uncertainty associated with 
the background pole. For the subthreshold state, the 
value of F^ strongly depends on the reduced width 7^. 
In fact, the product F^ x 7-^ does not depend strongly 
on the fit conditions [5]. Our values for F^ are mostly 
higher than the 90% confidence interval of 0.28-0.75 
eV from Ref. [6]. However, it is important to note 
that the >S'(0) value for the g.s. transition is not very 
sensitive to the fit parameters and we have adopted 
S{0) = 0.50 ± 0.05 keV-b. While this value seems to be 
at odds with that of Ref. [16], the 5(0) quoted there was 
for direct capture only and did not include the inter- 
ference effects of resonances, which are clearly important. 

Fits for the transition to the 6.18-MeV state have been 
performed, including El contributions from the 0.259- 
MeV resonance, the subthreshold state, and a back- 
ground pole. Figure 3 shows two fits with different values 
of F^ for the subthreshold state and different parameters 
for the background pole. Both fits give very similar (and 
large) values. We have adopted an intermediate value 
of S{0) = 0.04 ± 0.01 keV-b. Although the ^-factor for 
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TABLE III: Results of the i?-matrix fits for the g.s. transition 



a 

(im.) 


ps 

(e\ ) 


(K.e\ ) 


(^ev; 


5(0) 

I js-t; V u J 


X'/N 

(N — 91'! 


(keV b) 


10 

1 


3.5 


0.4 


1.3 


50 


0.52 


0.54 


4 


1.1 


1.3 


60 


0.53 


0.54 


o 




4.5 


1.7 


1.3 


60 


0.53 


0.54 




0.1 


5 


1.9 


1.3 


50 


0.47 


0.53 






5.5 


2.2 


1.3 


40 


0.46 


0.53 




0.01 


6 


2.6 


1.3 


40 


0.47 


0.54 




0.001 



J5 



0.001 




0.01 



1000 



E™(keV) 



FIG. 2: _R-matrix fit (a = 5 fm) of the present capture data 
to the ground state (closed circles). The corrected data from 
Sch87 are also shown (open squares) but are not used for the 
fit. 



this transition is rather uncertain, it has a small effect 
on the total 5-fac;tor. In Figure 3, we also show the 
Ml contribution calculated using the parameters given 
by Nelson et al. [17]. In contrast to their conclusion, we 
find this contribution to be negligible at low energies. 

The present extrapolated S-factor at zero energy is 
compared in Table V with previous results. In addition 
to the statistical uncertainty quoted there, we also esti- 
mate systematic uncertainties associated with the mea- 
surement of total charge (« 2.5%), photopeak efficiency 
(~ 3%), effective energy (« 3% at the lowest energies) 
and target composition (« 8%). If the systematic uncer- 
tainties are combined in quadrature, then they amount 
to an overall uncertainty of 9.4% at the lowest energies. 
Consequently, our result is S{0) = 1.64±0.07 (stat) ±0.15 
(sys) keV-b, which is 51% of the value of Sch87. Although 
our result for S{0) agrees with a rcanalysis of the Sch87 
data. [5] , the values for all 3 transitions arc quite different 
(and thus is the energy dependence of the S- factor). Note 
also that that our findings are consistent with another re- 
cent measurement of ^^N{p,^Y^O at low energies [8]. To 
further improve on the accuracy of S{0), new data are 
needed at lower energies, but also at energies above 0.5 
MeV. 

We have explored one of the consequences of this result 
by calculating the evolution of a star with a mass of 0.8 
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FIG. 3: J?-matrix fit (a = 5 fm) of the capture data to the 
6.18 MeV state (closed circles). The data from Sch87 arc also 
shown (open squares) but are not used for the fit. The solid 
(dotted) line is the fit with a background pole at 3.0 MoV (5 
MeV); the dashed curve represents an Ml contribution (see 
text). 



times the mass of the sun. The composition was appro- 
priate for the galactic halo, with an overall metalicity Z 
= 1.7 X 10~^ (or 1% of solar). With the present rate for 
"'^^N(p, 7)^^0 we find that the age at the main-sequence 
turnoff is 0.8 Gy older than that with the previous rate. 
The implication is that globular-cluster ages will have to 
be revised upwards. Further investigation of this issue 
and of other implications arising from this work are in 
progress. 



TABLE IV: Results of the i?-matrix fits for 6.18 MeV state 
for a = 5.5 fm 



ps 

(meV) 


^bkg 

(MeV) 


pbkg 

(MeV) 


pbkg 

(eV) 


5(0) 
(keV-b) 


{N = 19) 


5.0 


3 


1.5 


1.0 


0.04 


4.9 


0..-) 




1.-) 


0.1 


O.O:! 


■j.S 



TABLE V: Summary of 5(0) values (in units of keV-b) 



Transition 


Ref. [3] 


Ref. [5] 


Present 


(MeV) 




RC ^ 


1.55 ±0.34 


0.0811!:^;^ 


0.50 ±0.05 


KG ^ 6.18 


0.14 ±0.05 


o.o6l'>:l- 


0.04 ±0.01 


KG ^ 6.79 


1.41 ± 0.02 


1.63 ±0.17 


1.10 ±0.05 


Total 


3.20 ± 0.54 


1.77 ±0.20 


1.64 ± 0.07 
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